SUFFICIENT CONDITIONS FOR SAMPLING AND INTERPOLATION ON THE 

SPHERE 



JORDI MARZO AND BHARTI PRIDHNANI 

Abstract. We obtain sufficient conditions for arrays of points, Z = {Z(L)}l>i, on the unit 
sphere Z(L) C § d , to be Marcinkiewicz-Zygmund and interpolating arrays for spaces of spherical 
harmonics. The conditions are in terms of the mesh norm and the separation radius of Z(L). 



1. Introduction 

Let L p (E> d ) the Banach space of measurable functions defined in the unit sphere § d with 

15= / \f(z)\ p da(z) <oo, 



if 1 < p < oo, and 



oo = SUp |/(2)| < OO, 

ze§ d 



when p = 00. Here a stands for the Lebesgue surface measure in § d . 

For I > an integer, let Hi be the space of spherical harmonics of degree / in E> d i.e. the space 
of eigenfunctions of the Laplace-Beltrami operator on §> d 

A §d Y + l(l + d-l)Y = 0, Y eU h 

We denote by for L a nonnegative integer, the space of spherical harmonics of degree not 
exceeding L 

L 

I1 L = span (J Hi. 

1=0 

With respect to the inner product in L 2 (E> d ) the spaces Hi are orthogonal. We denote by hi and 
ttl the dimension of Hi and 11^, respectively. By Stirling's formula, txl ~ L d , when L — > 00. 
Let Yj 1 , . . . , Y^ 1 be an orthonormal basis of Hi- The reproducing kernel in 11^ is given by 

L hi 

K L (u,v) = = C d> ^i 1+A ' A) ((«^)), «,« G S", 

«=0 j=l 
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where (u,v) stands for the scalar product in d = 2 A + 2, C^lL - ^ 2 goes to a positive 

constant when L — > +oo and P^'^ are the Jacobi polynomials of degree L and index (a, 0), 
normalized so that 

To discretize the L p -norms in the space of spherical harmonics, we consider arrays of points 
on the sphere. More precisely, for any degree L we take tul points in E> d 

Z{L) = {z Lj G§ d :l<j< m L }, L > 0. 

This yields an array of points Z = {Z(L)} L>0 in 8 d . Denote d(u, v) = arccos(w, v) the geodesic 
distance between u, v G S d . 

Definition 1.1. Let Z = {Z(L)}l> be an array with > it l for all L. We call Z an L p - 
Marcinkiewicz-Zygmund array, denoted by LP-MZ, if there exists a constant C p > such that 
for all L > and Q G Hl, 

n-l m L p p ™l 

(i) ^ s < < d^fcEKM 

i/1 <P < oo, and 

sup < C sup 

wGS d j=l,—,rriL 

when p = oo. 

In other words, the L p -norm in E> d of a polynomial of degree L is comparable to the discrete 
version given by the weighted £ p -norm of its restriction to Z(L). For the unit circle, d — 1, the 
spherical harmonics are trigonometric polynomials. In this case, for = ttl = 2L + 1, J. 
Marcinkiewicz and A. Zygmund proved that the array of roots of unity form an L P -MZ array, 
[MZ37J, observe that C2 = 1. In higher dimensions the situation is more delicate. For = til 
there are no LP-MZ arrays when p 7^ 2 and the case p = 2 is open, see [|Mar071 . For big 
enough there are always L P -MZ arrays, see for example [MNWOOl FMllL 

Definition 1.2. Let Z = {Z(L)}l> be a triangular array with rriL < for all L. We say that 
Z is L p - interpolating if, for arrays {cL^L^o^j^mj of complex values such that 

sup — \°Lj\ p < 00, 

L>0 ^— ' 

J = l 

there exists a sequence of polynomials Ql G Ili uniformly bounded in LP such that 

Q(z Lj ) = c Lj , l<j< rrij. 

The concept of MZ and interpolating families is important in signal processing. Somehow 
these concepts are opposite in nature. MZ families are dense enough so that the L p -norm in E> d 
is comparable to a discrete version. On the other hand, interpolating families are sparse enough 
so that one can interpolate some given data. 
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Observe that Z is L 2 -MZ if and only if the normalized reproducing kernels of U L at the points 
Z(L) form a frame with frame bounds independent of L. Therefore, Z is L 2 -MZ when Z(L) is 
a set of sampling for II ^ with constants independent of L. Similarly, Z is L 2 -interpolating if and 
only if the normalized reproducing kernel of II^, fe^, at the points Z(L) form a Riesz sequence 
i.e. 
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(2) C-'J^Wlj 



2 < 



3=1 



da(z) < C |«ij| 2 , 
j'=i 



for any {a^j}^ with C > independent of L. When Z is both L 2 -interpolating and L 2 -MZ 
the normalized reproducing kernels at the points Z(L) form a Riesz basis, for more about these 
concepts see [Sei95l. 

We denote by d(u, v) = arccos(w, v) the geodesic distance between u, v e S d . 

Definition 1.3. An array Z = {Z(L)}l>o is uniformly separated if there is a positive number 
e > such that 

d(z Lj , z Lk ) > j——^, if J ^ k , 

for all L > 0. 

The left hand side inequality in (OQ) holds if and only if Z is a finite union of uniformly sep- 
arated arrays, also the L p version of the right hand side inequality in dU holds if and only if Z 
is uniformly separated, see HMar071 or HOCP111 for the general case of a compact Riemannian 
manifold. 

Definition 1.4. Let X be a subset ofE> d . The mesh norm ofX is 



p(X) = sup d(u, X) = sup inf d(u, z) 

ugS d u& d zeX 



The separation radius of X is 



8(X) = inf inf d(u,v). 

ueXvex\{u} 



The mesh norm of X C S d is therefore the maximal radius of a spherical cap which does not 
contain points from X, and the separation radius is the minimal distance between points in X. 

1.1. Main results. Our results in this paper are the following sufficient conditions for an array 
Z to be L P -MZ or //-interpolating. 

Theorem 1.5. Let 1 < p < oo and Z = {Z(L)}l> be an array in E> d such that for all L > 

S(Z(L)) > i 

where 9 > 2j\, j\ is the first zero of the Bessel function J\(t), and d = 2 A + 2. Then Z is 
L p -interpolating family. 
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Theorem 1.6. Let 1 < p < oo and Z = {Z(L)} L > be a uniformly separated array in S d such 
that for all L > 

(3) P{Z{L)) < ^, 

then Z is an LP -Marcinkiewicz-Zygmund array. 

We want to point out that there are other results about sufficient conditions for L P -MZ and 
interpolating families of points on compact manifolds, but such conditions do not provide precise 
constants, see HFM 1 1 1 IPMTOl IMNS W02l IOCP 1 2H . We observe that due to the result mentioned 
above about minimal L P -MZ (or maximal //-interpolating) arrays, an array with mi = 7Tl 
cannot satisfy the conditions of Theorems 1 1 . 6 1 1 .31 

When § 2 and for some particular arrays of points, there are some results about separation 
radius and mesh norm HRei90l ISW041 IDM 05 1 . The results we know are not very precise and we 
would just mention one to illustrate the use of our results. 

The set X = {x±, . . . , xn} C § 2 is said to be in s-extremal configuration if X maximizes the 
Riesz s-energy 

l<i<j<N 1 1 jl 

for subsets of N points on the sphere. In HKS98H it is assumed in order to get an estimate for the 
separation radius that te Voronoi cells around points in s-extremal configuration are all hexagons. 
Then it is obtained that 

The same way it seems reasonable to estimate the mesh norm by the value of the maximal radius 
of the hexagon getting 

>^0 1/Vi/2 

Therefore, if we take (kL) 2 points for degree L in order to assure we get an //-interpolating 
array we need k < 0.792, and to get an L P -MZ array we need k > 1.4. 

1.2. Outline of the paper. In Section 2 we prove Theorem 11.51 We use the classical approach 
by Ingham to obtain sufficient conditions for interpolation, [Ing36|. Ingham idea has been used 



in different context, HOU10I IKL05H . The main problem is the construction of appropiate pick 
functions, Lemma 12.11 which depend on an upper bound for the first eigenvalue of a spherical 
cap, HPin81llBCG83H . Estimates for the first eigenvalue are known also in general Riemannian 
manifols. 

In Section 3 we prove Theorem 11.61 Our approach follow the classical ideas of Beurling to 
study sampling sequences in Bernstein space, [Beu89 |. We define weak limits of an array, and 
relate uniqueness sets with the L°°-MZ property. Our result is consequence of a uniqueness result 
due also to Beurling. 

In what follows, when we write A < B, A > B or A ~ B , we mean that there are constants 
independent of L such that A < CB, A > CB or C\B < A < C 2 B, respectively. Also, the 
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value of the constants appearing during a proof may change but they will be still denoted with 
the same letter. 

2. Sufficient condition for interpolation 

In this section we prove Theorem 11.51 We adapt a nice idea of Ingham for Dirichlet se- 
ries, ]Ing36| . The idea rely on the construction of some pick functions with appropiate spec- 
tral properties. The existence of such functions can be established studying the first eigen- 
value/eigenfunction of the Laplace-Beltrami operator on the sphere from which a spherical cap 
has been removed, see HKL05llQUT0l Estimates on this first eigenvalue, HFH76llPmtniBCG83H . 
together with a result about perturbation of interpolating arrays, [MarOV , Lemma 4.11], provide 
the result. For the unit circle, d — 1, the condition in Theorem II .51 is 6 = n, and the proof is 
technically simpler. 

To prove Theorem 1 1.5 1 we use the functions given by the following lemma. 

Lemma 2.1. Given L and 9 as in Theorem U .51 there exists functions Fl such that 

(1) F L eL 2 ([-l,l]). 

(2) suppF L C [cos(0/L),l]. 

(3) [F L {{u, ■))}(£, j) < Oforall£> Land [F L {{u, ■))}{£, j) < 9 d/2 for all £ < L. 

(4) F L (1) ~ 7T L . 

where [F L ((u, ■))](£, j) stands for the Fourier coefficient J sd F L ((zLj, u))Yjf (u)da(u). 
Before establishing the existence of such functions we prove our main result. 



Proof. [Theorem [T31 Recall that the normalized reproducing kernel can be written as 



=0 fc=l 

Let F L {x) be given by Lemma [2T1 i.e. F L is a continuous function defined in — 1 < x < 1 such 
that 

suppF L c [cos(0/L),l], 
and F L (1) > n L ~ L d . Moreover, the Fourier coefficients of F L are negative for £ > L and 
uniformly bounded by CQ d l 2 for £ < L. Thus by using these estimates and Funk-Hecke formula 
we get 



J2 c Ljk L (z,z LJ ] 

3=1 



da(z) 



1 

L h 



1=0 k=l \j=l 

2 

m L 



da(z) 



u 1=0 k=l 
^ ^ +oo h e m L 

=0 k=l i,j=l 



~ 9 d / 2 7T L 
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mi +00 h? 

1 



ij=l fcO fc=l 

|cLj| , 



1 1 ^ 1 TP ^^\ L L 

1 1 x 1 rill) v~~\ 1 1 o v. \ ^ 1 1 2 



tj=i j=i i=i 

and Z is therefore L 2 -interpolating, because the other inequality in © follows directly from the 
separation. 

In order to prove the result for other p 7^ 2 we define, for 5 > 0, the per turbed array Z$ = 
{Z(L 1+S )} L , where L 1+s = [L(l + <5)]. It was proved in [Mar07„ Lemma 4.1 1] that if Z is L 2 - 
interpolating then Z_s is //-interpolating for all p E [1, 00]. Therefore, assume that Z satisfies 
the geometric separation condition 

i] := L min d(z u , z Lj ) > 9. 
Let S > be small enough so that t] > 9 + 59. We assume that L » 1 so that L5 > 1. Then 

Ld(z Ll+si , z Ll+d ) > -—77 > 0-^(1 + <5) > 

Thus, the perturbed array satisfies the same separation condition as Z, then is //-interpolating 
and Z = (Zs)s is //-interpolating for all 1 < p < 00. ■ 

Proof. [Lemma [2T]] Let C e /2L be the spherical cap of those x = (x\, . . . , Xd+i) G § d such that 
cos ^ < Xrf+i < 1. We denote by / the eigenf unction of the problem 

A§d/ + Vl/o = 

in C0/2L corresponding to the first eigenvalue Ao,l of the Laplace-Beltrami operator A§d. It is 
known, see [FH76J, that / belongs to the class of zonal functions, Lipschitzian, nonnegative, 
non identically zero and with support in [0, ^-]. 

We normalize H/0II2 — and define the zonal function 



(fo * fo)- 



L(L + d-l) 

with support in [cos f , 1]- 

Let {Yp}j : £ be the orthonormal basis in L 2 (E> d ) given by the spherical harmonics, then 

M(¥>) = EENv>Ej1*). 

e>o j=i 

where by Funk-Hecke 



[Fl((u,-)W,j)= F L ((u,v))Y l \v)da(v) = F L (i)Y l \u) 
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and 



here is the Gegenbauer polynomial of order a and degree t. 
On the other hand 

[Fl((u, ■))](*, J) = [ d ( 1 + L(L+ S d-1) ) (/o * fo)(M)y!(y)da(v) 



A S d(/ *fo){(u,v)): 



(fo*fo)((u,v))Yi(v)da(v)+ I 2(L+7-l) F/W0rW 



1 



mfYi{u) + L(L ^ - l (f * f )((u,v))A sd Y>(v)da(v) 
£{£ + d-l) 



/o(W(")- 2(L + d-l') / (fo*fo)(M)Yi(v)do(v) 



So we have proved that 



Note that the coefficients F L (£) < for all £ > L. Now we are going to prove that the 
coefficients fo(£) are bounded for £ < L. 



\fo(i)\ 





1 




-D/2 (1) 


< 


1 




- 1)/2 (1) 




AT 


~ r {d- 


-D/2 (l) 



1 



2 



/oW^ 1)/2 (t)(l-t 2 ) (d - 2)/2 rft 

e/2L \ 1/2 

L7f- 1)/2 (cos^) 2 sin d - 1 ^ 

\ 1/2 

L7^ 1)/2 fcos^ 2 sin d - 1 ^ < 



where we have used that C ( /~ 1)/2 (l) ~ £( d " 2 )/ 2 and cf ~ 1)/2 (i) < C^CV--*) (see [|Sze39l 
Section 7.32]). 

On the other hand, note that 

(/o * /o)(l) = (/o * fo)((N, N)) = ^ fg(cos9) sin*- 1 6d9 ~ tt l . 

./o 
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Thus, 

Fl(1) = (jo * /o)(l) + L(L+ 1 rf _ 1) A §d (/o * /o)(l) 
= (/o * /o)(l) + j - (L+ 1 rf _ 1) (/o * A sd / )(1) 

^ HL + d-l) ) (/o*/o)(l)-vr L (l- L(L ^_ i; 
So we need to find the smallest # so that the quantity 

1 ^ > 0. 

L(L + d-l) 

Equivalently, we need the smallest 9 so that Ao,l < L(L + d — 1). Using the upper bound from 
HBCG83H . we get that 

where ju-2)/2 is the first zero of the Bessel function J<d-2)/2- So taking 6 as in the hypothesis 
we have the result. ■ 

3. Sufficient condition for Sampling 

In this section we follow the classical approach used by Beurling to study sampling sequences 
in the space of bounded bandlimited functions, [Beu89]. First we identify the space of spherical 
harmonics composed with the exponential map as a subspace of the space of bounded bandlim- 
ited functions. Then we define the concept of weak limit of an array, and relate uniqueness sets 
with L°°-MZ arrays. Finally, we get Theorem II .6l by using a result of Beurling about uniqueness 
sets and a result about perturbation of MZ arrays, [Mar07. Lemma 4.9.]. 

Let x = (xi, . . . , x d ) G M. d and denote <fi(x) — (x\ H h x d ) x l 2 . The exponential map in § d 

is defined by 

sinc^fx) sva.4){x) ,, N \ ^ d 



exp(x) = xi , . . .,x d ,cos0(x) G 

V (p[x) <p[x) J 

Observe that exp(;z) is defined also for z G C d and is an entire function. 
Given Q L G II L we define the function 

(4) Q L (z) = Q L (exp(z/L)), zeC d , 

and the corresponding space IT L . Observe that 

sup \Ql(u)\ = sup \Ql{x)\. 

The following result shows that functions in YL L are entire in C d with Fourier- Laplace trans- 
form supported in the unit ball. 

Proposition 3.1. If Q G il^ then the Fourier transform, of Q has support in the unit ball ofM. d . 
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Proof. The reproducing kernel of Tl L centered at v G § d is, up to constants, the Jacobi polyno- 
mial u h-> P^ +X ' X \(u, v)). Let y G M. d be such that exp(y/L) = v. Consider the entire function 

C d 3z^ Pi 1+A ' A) (<exp(z/L),exp(y/L)>) . 
It is enough to see that for £ < L and for some constants C, N > (that may depend on L) 

\(exp{z/L),exp{y/L))f < C(l + \z\) N e^ z \ z G C d . 
For any C G C one has 2(3C) 2 = IC| 2 - ^C 2 - Therefore for any z = (z u . . . , z d ) G C d 

i=l i=l i=l 

and we get by the triangle inequality that 

l^i 2 + "- 




Finally 

\{exp(z/L),exp(y/L))\ = 



(z,y) sm<f)(z)/L . \y\ <j>(z) \y\ 

■ sin — + cos — — cos 



L\y\ 4>{z)/L L L L 



< 1 + 



( z ,y) ■ 1 2/1 

sin 



e l^)/ L l < C y>L (l + \z\)e 



L\y\ l 

and Q is the Fourier-Laplace transform of a distribution supported in the unit ball of M d . ■ 

Let B be the Bernstein space of entire functions in C d , bounded in M. d with Fourier transform 
supported in the unit ball of R d , endowed with the uniform norm. 

Given an array Z we send the points in Z(L) to M. d via the exponential map and define the 
corresponding family of weak limits. The Frechet distance between the closed sets A, B C W 1 
is given by 

[A, B] = mt{A C B + B(0, t), B C A + B(0, t)}. 

Definition 3.2. Let Z = {Z(L)}l>q be an array in E> d . We say that Ac! 11 is a weak limit of 
Z, denoted as A G W(Z), if there exist rotations pl G SO(d + 1), such that 

Lexp- 1 ( Pi Z(L))-A, 

where the above expression means that for any K C M d compact 

[(Lexp- 1 (p L Z(L)) n K) U <9fT, (A n K) U <9fT] ->■ 0, L ->■ oo. 

Proposition 3.3. Tfany A G W(Z) is a uniqueness set for B then Z is a L°°-MZ array. 

Proof. We argue by contradiction. Suppose that A G W(Z) is a uniqueness set for B but Z is 

not L°°-MZ. For any n G N there exists Q n G n Ln such that Q n {N) — ||<5n||oo = 1 and 

(5) -> SUp iQnfcOI- 

n j=l,...,mz,„ 
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From the sequence (Q n )n defined as in © it is possible to select a subsequence (see HNik75l 
3.3.6.]) converging uniformly on compact sets of C d to some function f E B with f(N) = 1. 
We denote this subsequence as before. For any A G A there exists a sequence Zi n j n G Z{L nk ) 
such that 

R d 3 w Lnk = L nk exp-\z Lnkjnk ) -»■ A, k oo. 
We denote this subsequence as before and we get 

|/(A)| < |/(A) - Q n (A)| + |Q n (A) - Q n (w Lnjn )\ + \Q n (w Lnjn )\. 

The first term on the right side clearly goes to zero. Also the last term goes to zero because of 
©. For the second term, using Bernstein's inequality we get 

|Qn(A) - Qn{z Lnjn )\ = |Qn(eX P (A/L n )) - Qn{z Lnjn )\ 

< L n rf(exp(A/L n ),2; Lnin )||Q n || 00 -> 0, n -> oo, 

because, see HBC731 p. 229], 

ci(exp y, exp f-) = ^ ~ ^ + o(L~ 1 ), x,?/GM d . 

LJ 1J LJ 

We get that / vanish in A, but f E B and therefore / = 0. ■ 

To prove our main result we use the following result about uniqueness due to Beurling, [Beu89 , 
p. 310]. 

Theorem 3.4 (Beurling). Let f be an entire function in <C d . Assume that 

hmsup =r<oo, / max log /(Oho < 00 • 

Assume f = on a discrete set A cM d such that 



r limsup inf \x — A| < — . 



TTjen / = 0. 



Proof. [Theorem \L6j Let Z be such that for all L big enough 

P(Z(L)) < i. 

Let A G and X L = Lexp' 1 (p L Z(L)) C R d with X L A. We want to see that A is a 

uniqueness set for B. 

For f E B and e > there exists A € > such that 

1/(01 <Aj 1+ ^, (eC d . 

Also | /(a;) | < M < oo for x E M. d , so we can apply Beurling's result (with r = 1) getting that 
any Acl d such that 

7T 

(6) limsup inf \x — A < — , 



is a uniqueness set for i3. 
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If 

sup ml \x — z\ < — , 

\x\<ttL z&X l 2 

for any L big enough, then we can deduce © which is equivalent to 

\T -It \ I 71 

sup mm |Lexp (cu) — z\ < —. 

u , eS dzeLexp~ 1 {p L Z(L)) ' 2 

But this follows from the condition on the mesh norm the property 

d{exp y , exp f-) = \^zA + o(L~ 1 ), x,y eR d . 

Therefore the condition on Z implies that it is an L°°-MZ array. In order to deduce the result 
for 1 < p < oo we define, for 5 > 0, the associated arrays Z$,Z-$ by Z$(L) = Z([(l + 
6)L]),Z_s(L)=Z (l(l-8 )L]). 

It was proved in [Mar07 , Lemma 4.9.] that if Z is an L°°-MZ array then Z$ is an L P -MZ array 
for all 1 < p < oo. 

Suppose that 

71 

V = supLp(Z(L)) < -, 

L>0 2 

and 5 > be such that 

r? < on. 

1 2 

For L big enough with Lb > 1 and any w e E> d we have 

L 7T 

Ld{u,Z Ll _ s ) < 7] < -, 

therefore Z^ 5 is L°°-MZ and (Z- S ) s = Z is L P -MZ. ■ 
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